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Abstract 

Stochastic Resonance in single voltage-dependent ion channels is investigated within a three state 
non-Markovian modeling of the ion channel conformational dynamics. In contrast to a two-state 
description one assumes the presence of an additional closed state for the ion channel which mimics 
the manifold of voltage- independent closed subconformations (inactivated "state"). The confor- 
mational transition into the open state occurs through a domain of voltage-dependent closed sub- 
conformations (closed "state"). At distinct variance with a standard two-state or also three-state 
Markovian approach, the inactivated state is characterized by a broad, non-exponential probability 
distribution of corresponding residence times. The linear response to a periodic voltage signal is de- 
termined for arbitrary distributions of the channel's recovery times. Analytical results are obtained 
for the spectral amplification of the applied signal and the corresponding signal-to-noise ratio. Al- 
ternatively, these results are also derived by use of a corresponding two-state non-Markovian theory 
which is based on driven integral renewal equations [I. Goychuk and P. Hanggi, Phys. Rev. E 69, 
021104 (2004)]. The non-Markovian features of stochastic resonance are studied for a power law 
distribution of the residence time-intervals in the inactivated state which exhibits a large variance. 
A comparison with the case of bi-exponentially distributed residence times possessing the same 
mean value, i.e. a simplest non-Markovian two-state description, is also presented. 

PACS numbers: 05.40.-a, 87.10. +e,87.15.Ya, 87.16.Uv 
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I. INTRODUCTION 

Stochastic resonance (SR) is by now a well-established phenomenon with wide 

spread applications in physics, chemistry, engineering sciences and the life sciences. It refers 
to the fact that in nonlinear stochastic systems an optimal level of applied or intrinsic noise 
can dramatically boost the response (or, more generally, the transport) to typically weak, 
time-dependent input signals. This fact plays a prominent role in biology with its abundance 
of a variety of threshold-like systems that are subjected to noise influences SR 
is often also characterized in terms of an underlying stochastic synchronization between an 
applied stimulus and an intrinsically stochastic dynamics . From the viewpoint of 

physical biology, the phenomenon of SR in biological sensory systems is commonly assumed 
to be rooted in the properties of excitable membranes. This being so, it ultimately can be 
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ained in terms of a driven stochastic dynamics of assemblies of ion channels 
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12t Il3j | . Although biologically relevant SR is generally a property of a cooperative 
coupling among ion channels Q, Q], the study of stochastic resonance in single ion 
channels carries merit on its own: (1) SR has not yet been convincingly demonstrated on 
the level of single molecules, with biological ion channels being such proper candidates Q; 
(2) ion channels can serve as suitable single-molecular sensors to be utilized in nanodevices. 

All these facts in turn have stimulated a vivid interest on this subject matter. Inspecting 
the recordings of ion current flowing through a single ion channel, as obtained from a typical 
patch clamp experiment bears close similarity with a simple, driven bistable dynam- 



xn clamp 

BBQ. 



ics |2|, [16( . The occurrence of SR in single ion channels thus seems very likely. In reality, 
the situation is, however, more complex, because (i) the dependence of the opening and clos- 
ing rates on the voltage and temperature are generally not Arrhenius-like 
a characteristic feature which plays a central role for the occurrence, or non-occurrence of 
the SR jsj], and (ii) the detailed studies of the statistics of the ion current switching events 
reveals that the probability distributions of the residence time intervals in different conduc- 
tance states are normally not single exponential. This implies that the observed dynamics of 
"on-off" conductance fluctuations is generically non-Markovian within a two-state stochas- 



tic description 2l|, [22]. In the simplest nontrivial case of a bi-exponential distribution of 



the residence time intervals dwelled in the nonconducting conformation, the e merg ing non- 
Markovian dynamics can be embedded into a three state Markovian dynamics with an 
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extra (visually indistinguishable) closed state. The presence of such a third state (or, more 
generally, a number of additional substates) can, however, be inferred from the bursting 
character of the observed dynamics when the ion channel dynamics after switching between 
its open and closed states for a number of times suddenly stops, and then revives again later 
after a notably longer time span has elapsed (as compared with the typical sojourn length 
in the open, or closed state). The presence of such a third state can be considered as a 
closed inactivated state. This "third" state in fact does not constitute a single state, but 
rather a manifold of many substates. As a consequence, the recovery of ion channel from 
its inactive state does not present an exponential rate process, but will be governed by a 
nonexponential distribution of corresponding residence times. This feature, i.e. the absence 
of a well defined time scale for recovery from inactivation can be observed in various ion 
channels (23 ] . 

Our main objective here is to investigate the basic features of SR occurring in single 
ion channels within the framework of a three state non-Markovian dynamics. In addition, 
we contrast this approach with the two-state non-Markovian theory of stochastic resonance 



developed in prior works 



24, 251. 



II. MODEL SETUP 



Following the reasoning put forward in Refs. [19, 20] we consider a discrete state model 
of ion channel gating, whose essentials are depicted with Fig. [T] It consists of three "states" : 
an open one (O), a closed one (C) and an inactivated state (/). The inactivated state / 
is assumed to correspond to the manifold of voltage-insensitive conformations of the ion 
channel protein and the opening C — > O and closing O — ► C transitions are associated with 
the motion of a voltage sensor and the opening/closing of the channel's gate. The stochastic 
transitions between the open and closed states are characterized by generally time- dependent 
opening and closing rates, k (t) and k c (t), respectively. These rates are assumed to depend 
on an externally applied, time-dependent voltag e sig nal. 
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we assume that in- 



In the spirit of the modeling put forward in [19J, I2C 
activation occurs from the closed state of the ion channel. In doing so, we are dealing with 
an archetype model of gating with opening, closing and inactivation dynamics Q]. Fur- 



thermore, unlike in the standard Markovian modeling 
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we assume that the transition 
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FIG. 1: Sketch of the three state model setup with generally time-dependent opening and closing 
rates and with a nonexponential residence time statistics in the inactivated state. 

from the inactivated, voltage-independent state to the closed state is not rate-limited, but 
rather is characterized by a broad distribution of rates. Put differently, we model the step 
from inactivation towards the closing state, cf. Fig. ^ by a non-exponential distribution 
t/>(t) of the residence time-intervals dwelled in the inactivated state. This distribution will 
be assumed to possess a finite average (r r ) := J °° Ttp(r)dT. It is voltage-independent but 
arbitrary otherwise. Furthermore, it is assumed that the channel's inactivation occurs from 
the closed state with a voltage-independent rate ki n . The voltage-independence of / <-> C 
transitions follows from the underlying character of conformational dynamics: Namely, one 
assumes that the I <-> C transitions occur in a direction being transverse to the direction 
of C „ O transitions ,20]. These latter transitions are ti g ht to the relocation of the gat- 
ing charge across the membrane, while the former transitions are not related to a charge 
redistribution; for further details we refer the reader to the discussion in Ref. j^oj. 



br several types of ion channels 
This circumstance in turn 



This chosen distribution of rates accounts for the fact that 
the recovery process does not have a well denned time scale 
implies that the gating process is non-Markovian within our three state description. The 
standard Markovian three-state description is recovered when t/>(t) = k r exp(—k r r), where 
k r := (tv) is the rate of the I — > C transition. 



III. NON-MARKOV THEORY OF THREE-STATE GATING 

The time evolution for the probabilities to occupy the open, closed or the inactivated 
state, p (t), p c (t) and pi(t), respectively, is governed by the following generalized master 
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equations: 



Po{t) = k {t)p c {t) - k c {t) Po {t) 



Pc(t) = -{k in + k {t))p c {t) + k c {t)p {t) + T(t- t') Pl (t')dt' , 

Pi{t) = k in p c {t) - [ T{t- t') Pl {t')dt' . (1) 

J to 

The Laplace-transformed kernel T{t) reads: 

f (s) = _^L, (2) 
l-V(a) 

with ^(s) denoting the Laplace transform of ip(r) ("tilde" denotes throughout the Laplace 
transform of a corresponding function). The solution depends also on the initial condition, 
chosen for example to read, p c {to) = 1, at initial time t = to- 
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, |32|, 33]. Alternatively, this very set of non-Markovian evolution eqs. can be obtained 
more directly as well: The terms not expected from naive grounds in Eq. $Q are the ones that 
involve memory. Let us assume that the ion channel is prepared in the state /, pi{0) = 1 
at to = and impose ki n — > 0, i.e. no returns are possible. Then, the leakage of probability 
Pi{t) is due to the transition into the state C, i.e. Pi{t) must equal (with such an absorbing 
boundary condition) the survival probability = J 4 °° ip{r)dT. With this leakage given 
by $(£) = — J ' T(t — t')&{t')dt' we readily find by use of the Laplace transform method the 
relation in Eq.(J2J). The description in Eq.Q portrays a driven (i.e. inhomogeneous in time) 
3-state , non-Markovian renewal process with the corresponding residence time densities 
(RTDs) HQ, 

readme; 

i/>o(t + r, t) = k c {t + r) exp (- jf" k c {t')dA , 

ip c {t + r,t) = [k in + k D {t + r)] exp (- [k in + k {t')]dt^j , (3) 

Note in particular that if) (t+T, t) and ip^t + r, t) depend not only on the length of residence 
time intervals r, but also on the entrance time point t. We refer the reader to Ref. [2^ 
for a detailed trajectory description of such driven renewal processes. These conditioned, 
nonstationary two-time RTDs ipj{r\t) := ipj{t + r,t) {j = o,c) are given as the negative time 



derivatives of the corresponding survival probabilities $j(r|t) := $>j(t + r,t); i.e. ipj{r\t) := 
-d$j(r\t)/dT. 

The averaged, time-dependent conductance g(t) of the considered ion channel reads 
(g(t)) = Po(t)9o + (Pc(t) + Pi(t))g c , where g a and g c are the conductances of the open and 
closed conformations, respectively. In the absence of a time- dependent signal, (g(t)) = 
(9) st = P S o9o + {p S c +pf)9c, where p]£ with a = o, c, I is the stationary solution of Eq. ((TJ) in 
the absence of driving. For the following we assume a periodic voltage signal given by 

V s (t) = Acos(nt). (4) 

The mean deviation of the channel conductance (Sg(t)) := (g(t)) — (g) s t from its stationary 
value thus reads at asymptotic times within a linear response approximation 



(6g{t)) = A\x{n)\coa{nt-<p(n)), (5) 

where (p(£l) is the corresponding phase shift l|, 13- 111 ©> x(^) denotes the linear response 
function in the frequency domain. The linear response result for the spectral amplification 



(SPA) of the signal rj(Q) is given by rj(fl) = \x(&)\' 



A. Linear response theory 

In order to evaluate the linear response function we expand the solution of Eq.Q as 
Pc(t) = YlkL-oo 1 "^) exp(—ik£lt) and p a (t) = J2T=-oc 1k(t) exp(— ikQt). The solution for 
Pi(t) follows by virtue of probability conservation, i.e., Pi(t) = 1 —p {t) —p c (t). In the limit 
t — > oo the solutions become tlm e-pe riodl c QQ with tae-independent ^efficients r t a,d 
q k which depend on the amplitude strength A and frequency Q. The asymptotic, nonlinear 
periodic solution thus reads 



PTit) = Yl r kexp(-tknt); r„ k = r* k , 

k=—oo 
oo 

P a o S (t) = Yl 9*exp(-iA;m); q. k = q* k . (6) 



The linear response function x(^) is encoded in the first Fourier term 



It reads, 



x(fi) = 2A^hm|, (7) 



where Ag := g — g c is the difference of ion channel conductances in the open and closed 
state, respectively. 

To start out, we assume the following small-signal expansion of the time-dependent re- 
action rates 

k 0>c (t) = is exp[-AG OtC (V(t))/k B T] 

« kM(l-p 0iC V s (t)), (8) 

where in absence of driving 

fc£> = u exp[-AG OtC {V )/k B T\, (9) 

with AG 0iC (Vo) denoting the corresponding static free energy barriers, Vo being a static 

voltage in the absence of signal, i.e., V(t) = V + V s (t), and /3 0jC = ^f^\ v=V() . 

The substitution of © into (JTJ) and taking (|1J) and (JHJ) into account yields for t — > — oo 
(this procedure is equivalent to taking to = and £ — ► oo) the following recurrence relations: 

[-ikSl + T(-ikQ,)](r k + q k ) = -k in r k + — -5 k0 , 

Vr) 

i-ikil + k^)q k + + g fc+ i) = fcf r, + fc«(r fc _a + r fc+1 ), (10) 

where 

In Eq. (110)1 . we used that lim s _>o r(s) = l/(r r ) which follows from the expansion ^>(s) = 
1 — s(r r ) + o(s), where o(s) stands for the terms such that lim s ^o — 0. From the first 
equation in (|TU|) we obtain 

q = 1 - (I + k in (r r })r (12) 

for A; = and 

q k = -(l+ — -g- — - )r fc (13) 



-ifcn + r(-ijfen). 

otherwise. Therefore, either g^, or can be determined from the second equation in Eq. 
()1U|) which thus uncouples into a recurrence relation for either q k , or r^. To the lowest order 
in A, q = q^ + 0(A 2 ) and r = + 0(A 2 ) where {go , r^} are the values when the 
perturbation is absent. A term linear in A is absent because a change from A — ► —A (which 



is equivalent to a phase shift by it) cannot result into a shift of steady state populations qo 
and ro- Moreover, q\ oc A and r\ oc A to leading order in A. Likewise, the expansions of q2 
and r 2 in A start out from A 2 . 

Then, by taking into account r = kcqo/ko + 0(A 2 ) in Eqs. (JTUJ), (fT2|) . (jTBjl we obtain 
after some algebraic manipulations 

r o = — } jfij- + 0{A 2 ) , 

1.(0) 

" ^.P^W) ^' (14) 



and 



91 = 5 A 7^X^T-fS L ~^> + 0(A2) ' (15) 

2 (r ) + (v) fc (o) _ in + 

iQ— ki n [l—ip(—iQ.)] 



where 



(To) = !/*?> (16) 

is the mean residence time in the open state. Furthermore, 

(r c ,) = (l + k in (r r ))/k^ (17) 

is the mean residence time within the set of closed states d = (C, I) in the absence of 
time-dependent driving. 

This latter quantity is defined as (r c /) := J °° Tip c i(T)dr, where ip c '(r) is the stationary 
RTD in the set of compound closed states, when no time-dependent signal acts, i.e. V s (t) = 0. 
This auxiliary quantity is obtained as follows: The channel is prepared in the closed state 
C at to = 0, i.e., p c {0) = 1 and the back transition O — > C is set to zero by imposing 
k c — > 0. Then, the solution of the first two equations in (£[]) in the absence of driving yields 
the stationary survival probability of the compound closed states as $ c '(i) = p c {t) + Pi{t) 
and ipc'if) follows as ip c '(j) = —d<& c '{r) / dr. Using this scheme we find, 



1 + 



" m s+f(s) 

and in virtue of ip c '(s) = 1 — s$ c /(s) 



*' (s) = zip fr - • (18) 



, (0) 

= , + - ' (19) 



s+r( s ) 



By use of Eq. (J2J) in Eq. (fHJj) this stationary distribution of the residence times can be 
recast as 

1.(0) 

Ms) = M = • (20) 

s + k® + k in (l-i!>(s)) 

Furthermore, (v) in Eq. (|T7|) follows from Eq. (fTSj) as (v) = $ c /(0). 
Let us introduce also the auxiliary function 

S(s) = d-toXi.-^W), (21) 
l - 4> (s)4> c i{s) 

where V> ( s ) = he /(s + k^) is the Laplace transform of the stationary RTD of the open 
time intervals. Then, Eq. (JT5j) can be rewritten in a more compact form as 



^l A /;:<\ Si ~ i o ) +o w- < 22 » 

2 ( r o) + (v) 

Since k in (T r ) does not depend on voltage V, /3 Q / := dl ^ c/ ^ | y=v = — dl ^y \v=v = Po and 

re result of the phenomenological two-state theory of 



the result in Eq. ()22j) coincides with t 



non-Markovian SR put forward in 
the linear response result 



24|,|2i|. Upon combining with Eq. ©, Eq. (J22) yields 



m = (& -A)^g(-«) , (23) 



(r ) + (p) 




in Ref. I25J. 





The just outlined procedure can as well be extended into the nonlinear response regime to 
obtain nonlinear response functions of required order in the signal amplitude. This, however, 
is beyond the scope of this work. This central result in Eq. (j23|) is alternatively derived in 
Appendix A by making use of the two-state non-Markovian theory detailed in Ref. j^ . 

B. Explicit results in terms of thermodynamic free energies 

We next introduce formally the effective free energy bias e(T), i.e., 

W - ■ (-) 



(v) V. k B T 

In accordance with the relations (Jlbl IT7j) we obtain 



e(T) = AG o (V ) - AG C (V ) + k B T\n(l + k in (r r )) . (25) 
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Using that AG o (V ) = G* — G c and AG C (V ) = G* — G , where G* is the free energy of the 
transition state, and G afi = H 0fi — TS 0fi is the free energy of the open (closed) conformation, 
Eq.()25|) can be recast as 



e(T) = AH- TAS + k B T\n(l + k in (r r )) , 



(26) 



where AH := H a — H c is the difference of thermodynamic enthalpies of the open and closed 
conformations (or of the internal energies, if no volume change of the macromolecule occurs 
at the conformational transition) and AS" := S a — S c is the corresponding entropy difference. 

Assuming that the free energy barriers in Eq. (jHj) have linear dependence of voltage, 
AG {V) = AG - q(l - 8)V and AG C (V) = AG C + qd V, where q is the gating charge 



and < 5 < 1 is a constant 



14. 



we find that (3 C — (3 Q = —q/{k B T). For the spectral 



amplification of the conductance response rj = |x(fi)| 2 we obtain then from Eq. (|2"H|) our 
first main result 

(Ag)V v\T) |£(-^)| 2 



16(fc B T) 2 cosh 4 



2k B T 



9? 



(27) 



with e(T) given in Eq. (PEJ, while u(T) = (To)- 1 + (r^) -1 denotes the sum of the closing 
and opening rates. Q{s) in Eqs. (|2T|) . (f2T|) reads explicitly 



s + + 



sk 



(0) 



(28) 



s+k in (l-il)(s)) 

where kf\kf\k in are the unperturbed rates and ip(r) is the residence time distribution. 
The result in Ea. (l27|) coincides, apart from a model specific constant, with the expression 
given in Ref. |2J], cf. Eq.(25) therein. Within the linear response approximation the 
corresponding signal-to-noise ratio (SNR) at the angular driving frequency Q is obtained 
from Rsn(Q,T) := ttA 2 !]/ Sn(uj = Q), where Sn(u) denotes the spectral power density of 
conductance fluctuations in the absence of signal V s (t). By use of the Stratonovich formula 
for the autocorrelation function of the alternating renewal process |35f and the Wiener- 
Khinchin theorem one finds for Sjy(u) 25l l35| 



S n (uj) 



2(Agf 



1 



(r ) + (v) ^ : 



-Re 



Q(iuj) 



(29) 
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Therefore, the SNR equals the result derived in Refs. J24j, |25j, namely, 

nA 2 q 2 v(T) 



Rsn(^, T) 



8(fcsT) 2 cosh 2 



2k B T 



(30) 



where 



N{tt) 



(31) 



provides the specific function which accounts for manifest non-Markovian effects. In the 



24, 



251 ]. where 



low frequency limit fl — > 0, N(Q) approaches the limit iV(0) = 2/(C 2 + C 2 ) 
C and C c i are the coefficients of variation of the RTDs of the open and compound closed 
states, respectively. 

By use of the expansion ip(s) = 1 — s(r r ) + o(s) in Eq. (|2~%j). and Eq. (JT7J) we find in 
the adiabatic limit Q — > that the spectral power amplification acquires the universal form, 
reading 

{Agfq 2 1 



lQ(k B T) 2 cosh 4 



2k B T 



(32) 



This result holds in the presence of asymmetry with nonvanishing e(T). In such a case, 
the SPA f!32|) can exhibit a sharp stochastic resonance at the physiologically relevant tem- 
peratures when e(T) changes sign, when the probabilities of the channel to be open, or to 
stay close become equal. This corresponds to the opening threshold for the detailed below 
case of an ion channel sensitive to cold (or, vice versa, sensitive to heat, if the open state 
is preferred from an entropic viewpoint). The cold-sensitive, or heat-sensitive ion channels 
36j | present appropriate candidates to reveal this entropic-SR effect. This feature at very 
small driving frequencies, which mimics a Markovian behavior, has thus nothing to do with 
non-Markovian properties; it is solely due to an entropic asymmetry. The non-Markovian 
effects emerge, however, at small, but finite frequencies in the corresponding SPA-curves, see 
in Ref.j^]. Most prominent non-Markovian, l ong- time memory effects distinctly suppress, 
however, the SNR in the low-frequency domain [2A 

7rA 2 q 2 



R SN (^^0,T) 



251 ] . For f2 — » 0, we obtain (with C Q = 1) 
u(T) 2 



8 (Wcosh 2 



2k B T 



1 + c 2 



2 • 



(33) 



with C c i > 1. Thus, the Rsn(Q = 0,T) is fully suppressed, making the detection of low- 
frequency signals barely possible, e.g., for the power law distribution iP(t) oc l/r 2+a with 
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< a < 1 considered below. This suppression occurs due to noise feature in the 

spectral power Sjv(^) for small frequencies. Let us illustrate now these general considerations 
with two particular models. 



IV. ION CHANNEL GATING: BI-EXPONENTIAL DISTRIBUTION VERSUS A 
POWER LAW 

The current three state model provides a suitable framework to clarify the role of a power- 
law distributed residence times in the inactivated and closed states as compared with the 
simplest two-state non-Markovian situation (with respect to the observable dynamics) of a 
bi-exponential distribution, that, likewise, can be embedded into a three-state Markovian 
description. 



A. The case of a bi-exponential distribution 

We start our driven channel gating investigation with the simplest case of an exponential 
residence time distribution of the transition I — > C, i.e., ip{r) = k r exp(—k r r). Then, Eq. 
flUJ) yields 

(s + ko )(s + k r ) + k in s 
The inversion of Eq. (JMj) yields a bi-exponential probability density, i.e., 



Us) = , . , ; , . • ( 34 ) 



■0c'(Y) = ciAi exp(-Air) + c 2 A 2 exp(-A 2 r) (35) 



with the rate coefficients Ai >2 = \{k^ + k in + k r ± y (k^ — k in — k r ) 2 + Ak o ^k in ) and cor- 
responding weight factors c i)2 = \ {\ ± (k ^ — k in — k r )/ \J (k^ — k in — k r ) 2 + Ak^ki^j . 
The mean residence time corresponding to this RTD (|33j) is 

(V) = 4- (l + k -P-) (36) 



ko 



(°) V k 



Such a non-exponential, i.e. bi-exponential RTD can actually be very broad as characterized 
by the corresponding coefficient of variation C := a/ (t 2 ) — (t) 2 / (t), yielding from 



2k o 0) k m , . 

^ = ^ 1+ ffc - +fcV - (37) 
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Indeed, for k r <C k in <C ko we deduce from Eq. (JHTJ) that « \j2k { o ] /k in > 1. This 
implies that the distribution (J35)) has a small, but very broad long-time tail, which in turn re- 
sults in a large variance of the residence times. This finding carries important consequences: 



As shown in Refs. 24, |2£j the signal-to-noise ratio Rsn is then strongly suppressed in the 
low-frequency limit Q — ► by the factor 1/N(0) = |(C 2 +C^), where C is the coefficient of 
variation for ?P (t), C q = 1 in the present case. This presents a first manifest non-Markovian 
effect which is present already within this simplest non-Markovian setting. We also note 
that the auxiliary function Q{s) can be recast as 

c( a \ — g ( g ~j~ kin + k r ) , , 

&W ~ [s + fii(T)][s + (j, 2 (T)Y (38) 

where 

/i lj2 (T) = l - (k® + fef + fc in + k r ± ^ (fcf + kf ] - k m - A; r ) 2 + Ak^K^j (39) 

are the decay rates of the conductance time-correlations. By use of Eqs. flHH|) . (j3Tj) the 
results in Eqs. (J2*T|) for the SPA, and (|3TIj) for the SNR assume explicitly the form: 

(Ag)y v\t) tf + {k m + k r y 



2k B T 



and 



R fOTl vrA 2 g 2 „(T) ^ 2 + {k in + k r f 

K SN (^,ij = — — — — — , (41) 



8 (^T) 2 cosh 2 [jgl] fi 2 + (A; . n + fcr)2 + k (0) k 



in 



respectively. The remaining parameters are: u(T) = (r ) 1 + (r c /) 1 = k^ + k^k r /(k r + fcj n ) 
with fc$ in Eq. © and e(T) = AG o (V ) - AG C (V ) + k B T \n(l + k in /k r ) with exponential 
Arrhenius rates k in = z/ exp(— AG in /k B T) and k r = v exp[— AG r /ksT] |37| . 

The results in Eq . (|4U|) and (|41j) constitute central results for the SR occurring in a three 
state Markovian model of gating in ion channels possessing an inactivation from the closed 
state. At the same time, these results correspond to a simplest non-Markovian two-state 
model of the observable dynamics of conductance fluctuations. 

We performed numerical calculations for the set of test parameters given in the Table 
which is chosen to mimic the experimental temperature dependence of the cold-sensitive ion 
channels, see in Ref. 
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TABLE I: Free energy barriers, AG a = AH a - TAS a 



a Enthalpy part AH a , kJ/mol a Entropy part AS a , ks 

c 175 52 

o 15 -10 

in 15 -15 

r 15 -20 



"Table depicts the corresponding parameters for the free energy barriers that enter the related rate coeffi- 
cients. Here, AH a must be divided with the Avogadro number Na to obtain the corresponding value per sin- 
gle molecule. For example, AH C — 175 kJ/mol thus corresponds (we use the Boltzmann constant fee and not 
the gas constant R = UbNa in the rate expressions, respectively) to a value, AH C w 2.91 TO -19 J « 1.81 eV, 
etc. 

The temperature-dependence of the corresponding Markovian transition rates are given 
in Fig. 2(a). For this set of parameters, the coefficient of variation of the closed residence 
time (compound state) within the three-state Markovian description is C c i ~ 17.14 and 
the low-frequency SNR is suppressed by the factor of l/iV(0) ~ 147.43 as compared with 
the corresponding two-state Markovian case with the same rates k^ and k^ := {t c ')~ x - 
Its temperature-dependence is depicted in Fig. 2(b). The frequency prefactor z/ has been 
taken to be v Q = 6.11 • 10 12 s _1 . This value corresponds approximately to a gas-phase 
value of ksT/h at T = 20°C It clearly overestimates any effects that relate to friction and 
other details determining the transmission coefficient k in condensed phases. The correct 
estimation of this prefactor would require a more elaborate theory of the Kramers type 
l^l rather than Absolute Rate Theory used here. Nevertheless, this ambiguity does not 
play role if to assume that the frequency prefactor uq is one and the same for all the rate 
parameters. This is so because a different value of z/ would result in one and the same 
(unknown) systematic entropy correction for all metastable states. The entropy differences 
between the metastable states remain unchanged. 

Our results for this choice of parameters are depicted in Fig. 3(a) for the spectral ampli- 
fication rj, and Fig. 3(b) for the corresponding signal-to noise ratio Rsn- In particular, Fig. 
3(a) convincingly demonstrates the possibility of Stochastic Resonance within the range of 
physiological temperatures when e(T) w 0, in accordance with Eq. (J3~2~j) . The SNR in Fig. 
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FIG. 2: Temperature-dependence of the rates used for the ion channel gating dynamics, absent the 
driving. Part (a) for the considered three state Markovian model; part (b) for the corresponding 
two-state non-Markov model with the effective opening rate defined as kjS := (t c /} _1 , where (r c /) is 
the mean residence time. The rate denotes again the transition from the open state O towards 
the closed state C in the absence of time-dependent driving. 

3(b) also does exhibit the SR-like increase with increasing temperature. It does not, how- 
ever, reach a resonance peak-behavior in the corresponding temperature range. This peak 
behavior would formally be assumed (not depicted with Fig. 3(b)) only at physiologically 
unrealistic high temperatures of about 600 °C. This resonance behavior is formally contained 
in Eq. ()33j) (notice the presence of v(T) in the numerator there, which invalidates the cri- 
terion e(T) 0). The behavior at extreme high frequencies Q 3> k in , k r , k^ (cf. Eq. (1411) ) 
assumes qualitatively the same behavior as depicted in Fig. 3(b); the only difference being 
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temperature [°C] 

FIG. 3: Markovian three-state model of ion channel gating: (a) spectral power amplification i] (in 
units of (Agq) 2 /(J/mol) 2 ) and (b) signal to noise ratio Rsn (in units of ir(qA) 2 /(J/mol) 2 ■ s _1 ) 
versus temperature (in °C) at different angular frequencies Q of the harmonic signal. In Fig. (a), the 
full line corresponds to the adiabatic limit Q — > 0; the dash-dotted line corresponds to fi = 10 s _1 
and the dotted line corresponds to Q = 100 s _1 . The corresponding lines for = 0.1 s _1 and 
{} = 1 s _1 cannot be resolved from the zero-frequency limit. In part (b), all the lines for the 
above mentioned four different frequencies merge with the zero-frequency result (solid line). Note 
that for the two-state Markovian counterpart of the considered three-state model with the rates 
depicted in Fig. 2(b), all the lines for the different frequencies merge with the zero- frequency result 
in part (a). Likewise, in part (b), the depicted line becomes then be multiplied with the factor 
l/iV(0) ~ 147.43 to yield the corresponding Markovian two-state result. 
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that the SNR becomes increased by the factor (C 2 + 1)/2. Surprisingly (in view of the rather 
large enthalpic barriers), this does not exclude the appearance of Stochastic Resonance in 
the SNR-behavior at physiological temperatures at some intermediately high frequency of 
the signal. For the studied model system this occurs, for example, for Q = 10 4 s _1 (this 
feature is not depicted because the accompanying spectral amplification rj(T) in Fig. 3(a) 
emerges to be extremely small). 



B. The case with a power law distribution 

To model the gating dynamics more realistically we next consider a description with 
a nontrivial nonexponential residence time distribution in the inactivated regime ip(r) in 
Eq.([28|). In particular, we use a probability density i/j(t) with the following characteristic 
function [29| 

^(s) = i / — \ — 7 V > ( 42 ) 

rK ' l + s(T r )g a (sr d ) v ' 



where 



tanh(z a/2 ) 

9a(z) = z { a/2 ] - (43) 



The RTD (1421) has been obtained in Ref. 



•3 



as solution of a conformational diffusion 
model accounting for an anomalous subdiffusion over energetically quasi-degenerate sub- 
states within the given conformation, i.e. / in the present work. This probability density 
possesses three parameters: the mean residence time (r r ) in this inactivated conformation, 
the conformational diffusion time Td and the index of subdiffusion a, < a < 1. 

In spite of possessing three independent parameters only, the RTD (J4"2~j) is capable to 
display a rich behavior. The case of single-exponential distribution with the rate parameter 
k r := (T r ) _1 is rendered for Td = 0. The case a = 1 corresponds to normal diffusion. For 
Td *C (r r ), the intra-conformational diffusion effects are not essential. However, for r d ^> (r r ) 
and in the range (r,.) 2 /^ <Cr<r^ the corresponding RTD assumes a negative power law, 
i.e., ip(r) oc r~ 3 / 2 , which ends up in an exponential tail for r > For a < 1, the 

distribution (}4"2*|) has an infinite variance, since (r 2 ) = oo for any Td ^ and, depending on 
a subtle interplay of parameters, it can display up to three different power laws. These are: 
ip(r) oc r~ Q / 2 initially, ip(r) oc t~( 2 ~ q / 2 ' ) intermediately, and V ; ( r ) oc r~( 2+Q! ) asymptotically 
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3- 



291 ] . With the discussed choice of tp(r), the characteristic function of the RTD in the 
compound C state reads 

ko°\sg a (sT d ) + k r ) 
s + k ( o y )[sg a (sT d ) + kr] + k in sg a (sr d ) 
It possesses the same average residence time (|36J) and the coefficient of variation is 



oo, a < 1; 

Cc = 4 / on ( 45 ) 

P^ + W, a = l. 

Note that for a < 1, C c / = oo. This implies that the detection of low-frequency signals SI — ^ 

is fully suppressed as Rsn{Q) — > at Q — > 0,13- A bistable stochastic element with such 

properties can thus be used as a high-pass filter for the signal transduction. In this case, the 

first leading terms of the small- s expansion of ip c '(s) read ip c '(s) ~ l - ( T <^) s ^'^~w a — (sr d ) 1+Q . 

This corresponds asymptotically to a distribution, Vv(t) oc r~( 2+a \ similar to the Pareto 

law behavior considered in Refs. 2^, 13 • 

With (J4*2~j) in Eqs. (|2*5ji. (J2*T|) and (jSOJ) one can evaluate the spectral amplification of signal 
and the signal-to-noise ratio. We did this for the above set of thermodynamical parameters 
(see in Table U} entering rates which is supplemented with the following parameters of the 
RTD in Eq. (|42j) : a = 0.25, r d = 0.1s. As it can be seen in Fig. 4(a), the non-Markovian 
effects resolve the lines at the frequencies Q — 0.1 s _1 and Q — 1 s^ 1 which merge with 
the zero-frequency line in Fig. 3(a). Moreover, the different lines which merge in Fig. 3(b) 
become also resolved, see in Fig. 4(b), where upon increasing the angular frequency Q the 
SNR Rsn{Q,T) grows. Namely, at zero frequency Rsn(Q — ► 0, T) = and for Q rs 100 s _1 
Rsn{&, T) reaches approximately (from below) its low-frequency three state Markovian limit 
in Fig. 3(b). With the further increase of Q SNR will, however, grow further approaching 
asymptotically the two state Markovian limit, where the non-Markovian form factor assumes 
unity, i.e. N(Q) = 1. It is worth noticing that for the two lower curves in Fig. 4(b) the 
SNR seemingly saturates with increasing temperature. In fact, however, the SR-behavior 
exhibits a rather broad maximum. The occurrence of this maximum is quite surprising (see 
the discussion at the end of the subsection IV. A) and is due to non-Markovian effects. 

Moreover, at variance with the SNR behavior (where the SNR increases with increasing 
angular frequency of signal) the corresponding SPA diminishes with increasing angular fre- 
quency. Therefore, we do not observe an intermediate frequency regime which would prove 
optimal for the detection of non-Markovian SR. 
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FIG. 4: The case with a power law distribution of residence times for three-state ion channel 
gating, cf. Fig 1. Instead of a single-exponential RTD ip(r), the residence time distribution in the 
inactivated state is given now by Eq. (|42j) with a = 0.25 and = 0.1 s. The other parameters 
remain the same as in Fig. 3. The manifest non-Mar kovian effects result in a resolution of the 
different lines that merged in Fig. 3. Note in part (b) that the signal to noise ratio is fully 
suppressed towards zero for Q — > (the result merges with the horizontal axis). The adiabatic 
limit of the corresponding Markovian three-state model is approximately assumed at $7 = 100 s _1 ; 
cf. the dotted line in (b) which compares with the solid line in Fig. 3(b). 

V. RESUME 

With this work we have investigated Stochastic Resonance in a three state non-Markovian 
model of ion channel gating. We note that our scheme of a three-state non-Markovian mod- 
eling is distinctly different from a similar one, recently applied to an excitable neuronal 
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dynamics |38|, |39|. The latter model assumes a three state system which cycles unidirection- 
ally between three states, 1 (silent), 2 (excited) and 3 (refractory), i.e., 1 — > 2 — > 3 — > 1. 
This unidirectional cycling corresponds to processes that are very far from the thermal equi- 
librium and which require a continuous supply of free energy. In contrast, our modeling is 
compatible with the thermal equilibrium. The ion channel gating is commonly assumed to 
be a thermal equilibrium process This is in contrast to the situation with ion pumps 
and neuronal systems which do require a free energy supply for proper functioning. As it 
has been shown in this work, our three state description is compatible with the phenomeno- 
logical two-state theory of non-Markovian Stochastic Resonance put forward in Ref. J^j], 
while the non-equilibrium, three-state model of Ref. 39] is not within this class of system 

n n 

behaviors |25|]. Nevertheless, the result of Ref. [39] for the spectral power amplification 



can be reproduced from our Eq. (jAl|) by specifying the RTD of the compound silent state 
(i.e. the silent state plus the refractory state) as a time-convolution of two corresponding 



residence time probability densities and expanding it like in Eq. (|A2j) . In other words, the 
model of Ref. Q ean be reeovered as a special ease within our general noneqnilibrinm 



approach put forward in Ref. 25], cf. section IV therein. 

In conclusion, the non-Markovian SR effects such as the suppression of the signal-to-noise 
ratio for low-frequency signals and the resolution of different signal frequencies in the spectral 
amplification of signal in the presence of possibly large (entropic) asymmetries can be nicely 
modeled already within a three state Markovian model. Such a three state Markovian 
model yields the simplest non-Markovian model after the corresponding projection onto the 
subspace of two observable states. The signal-to-noise ratio remains, however, finite in the 
limit of the zero-frequency signal. Its complete suppression requires an infinite variance of 
the residence time intervals in the inactivated state. For such a manifest non-Markovian SR- 
behavior, weak external oscillating signals with an intermediate frequency should be used 
in order to detect SR experimentally. 

Our present results provide a theoretical proof for the occurrence of Stochastic Resonance 
in single biomolecules at physiological temperatures. This being so, our findings can guide 
the experimentalists to choose both, the appropriate molecules for doing experiments and 
to identify the corresponding experimental parameter regime which in turn will reveal the 
SR-phenomenon in a single ion channel. Presently there exist only very few experimental, 
ion channel-based SR-works: In Ref. []| attempts to identify SR in a single ion channel 
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have been made; the occurrence of SR therein is, however, not convincing. In Ref. |6| one 
investigates SR and demonstrates SR. The phenomenon has been studied, however, only on 
the level of a small number of dynamically self-assembled alamethicin ion channels. 

Our findings show that ideal candidates to observe experimentally SR on the level of single 
molecules are cold (or heat) sensitive ion channels j36|. From this viewpoint, we believe that 
a repetition of the experiments in Ref. []| by resorting to such channels would indeed become 
successful for observing SR in a single ion channel operating within its physiological regime. 
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APPENDIX A: DERIVATION OF LINEAR RESPONSE FUNCTION WITHIN 
NON-MARKOVIAN TWO-STATE THEORY 



In this Appendix the result in Eq 

n 

theory of Ref. [25|| . Namely, in Ref. 



33J) is re-derived from the non-Markovian two-state 
25| it has been shown that the linear response function 
of an arbitrary two state renewal process characterized by the conditional RTDs ^ J=12 (r|t) 
and the amplitude of conductance fluctuations Ag to a periodic signal in Eq. (@J is given in 
the limit A — > by 



1 



2iAg 

'~An(n) + (r 2 ) 



i-^i(-in)^2(-in) 



(Al) 



In Eq. (jAip . if)j(s) denote the Laplace-transforms of stationary the RTDs in the absence 
of driving and ipj (s) are the Laplace-transforms of the corresponding contributions in the 
expansion of the conditional, driven RTDs, i.e., 

oo 

'UAt)= E ^ n) (r)exp(-innt). (A2) 
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The quantities V>j^(s) must be evaluated from an underlying multi-state, or continuous state 
dynamics to first order in A. Whenever ^^(s) satisfy the relation 

j>?\-iO) = -~/^4[l - ^(-m, (A3) 

then the result of phenomeno logical theory in Eq. (J2*3|) is recovered from Eq. (jAljl . 
For the exponential form of distribution of open residence times ip {r\t) = ^(t^) = 
k c (t + r) exp[— f* +T fc c (r')c?r'], the relation (jA3|) is fulfilled. We identify here the state "1" 
with the compound closed state "(C, I)" and the state "2" with the open state "O"; note 
also that (3\ = j3 , $2 = @ c - This relation is valid as well for a multi-exponential distribu- 
tion which assumes a scaling relation among the rate parameters which is not modified by 
driving, i.e. a form- invariant RTD, cf. j^. Below we demonstrate that the relation (|A3|) 
is also valid in the present case which obviously does not belong to the latter universality 
class. Nevertheless, the application of the phenomenological linear response theory to the 



considered non-Markovian processes 2M ls justified. 

In order to obtain the conditional survival probability of the compound closed state 
<3>i(t|£) := + r, t) one needs to solve the system of equations 

^Pc(r\t) = -[k in + k {t + r)} Pc {r\t) + J T{r - r') Pl {T'\t)dr' , 

^pPi(r\t) = k mPc (r\t) - jT r(r - T') Pl {r'\t)dr' , (A4) 

with the initial conditions p c (0 |t) = 1 and P i(0\t) = 0. Then, $i(r|t) = Pc (j\i) + 
P i(r\t). Asymptotically, the considered periodically driven renewal process becomes cyclic- 
stationary. For such a cyclic-stationary process, $i(r|t) must be invariant under time shifts 
t —>■ t + nT, n = ±1, ±2, ... with period T = 2tc/Q; it thus can be expanded into Fourier 
series like in (jA2|) . With k Q (t) being a periodic function of time, we are seeking a solution of 
Eq. (jH in the form, Pa (r\t) = EZ-ooP^\ T ) exp(-mfit), pi~ n) (r) = [ P ( a\r)}*, a = c,I. 
Invoking additionally a small-signal expansion (JHJ) this yields an infinite system of coupled 
integro-differential equations for (r) which upon the use of the Laplace-transform method 
results in a recurrent-difference relation for p^\s), i.e., 

8 + hn + ~ ^S^] P <?\s) + k^[ P ^\ S - in) + P { r 1] (s + «])] = S n , , (A5) 

s + r(sj-i 

and a relation expressing Pj(s) throughputs), reading, 

Pf\s) = ^—/:\s\ (A6) 
s + Y[s) 
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where oc A is given in Eq. ((lip . Note that p« (s) are nonlinear functions of the driving 
amplitude A. The equation ()A5|) can be solved perturbatively by using the corresponding 
expansions of pa(s) in A. The expansion of pa\s) in A starts from the driving independent 
terms and obviously (from the symmetry reasons, cf. A — > —A) does not contain linear 
contribution in the amplitude strength A. Moreover, p^^^s) oc A, to the lowest order in 
A. Hence, p^ 2 \s) = 0(A 2 ) because the response at the second harmonic driving frequency 
cannot be linear. The use of standard perturbation theory then yields to the lowest order 
in A, 

ms) = , + k J°> ^ +0(A2) - (A7) 

S -f K in -f /Co fi+f (s) 

Eq. ()A7|) determines in combination with Eq. ()A6|) for n = the stationary (i.e., in the 
absence of driving) survival function of the compound closed state in Eq. (|18|) and the 
corresponding RTD (JT9J) . respectively. 

Eq.flHl) yields in virtue of p c 0) (0) = l/ki° } + 0{A 2 ) and Eq. PJ the result 



Together with Eq. ()A6|) for n — 1 this gives with (3 Q = Pi 

$W(-iO) = ^x^f^-zQ). (A10) 

Upon using the general relation, ifia(s) = 8 Uj o — s$a (s), Eq. ()A10|) thus yields the desired 
relation in (1 A3)) . 
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